Effect of Electron Correlation on the Bragg Reflection 
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We study the effect of correlation on the Bragg reflection in the 3D electron gas, the ID Luttinger 
liquid, and the ID Hubbard model in an alternating periodic potential at half-filling. In the last 
system, we suggest a Luttinger-liquid-type quasi-metallic state in the crossover region from the 
band insulator to the Mott insulator. We explain the appearance of this state in terms of the 
incompatibility of the Bragg reflection with the concept of Luttinger liquids. 
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One of the fundamental issues in solid state physics 
is to elucidate to what extent the concept of "electronic 
band structure" is relevant in a strongly-correlated sys- 
tem. In the free-electron gas, the band gap is formed by 
the quantum interference between an incident electron 
plane wave, e lk r , with a reflected one, e I ( k + K ) r j where 
K is a reciprocal-lattice vector related to periodicity of V 
an electron-lattice potential. In this sense, the quantum 
coherence leading to the Bragg reflection is a key to the 
band structure and thus we consider it quite important 
to study the correlation effect on the Bragg reflection. 

The study is conceptually simple in the Fermi liq- 
uid as represented by the three-dimensional electron gas 
(3DEG) at metallic densities. The point is that we should 
grasp the band-gap formation in terms of the Bragg re- 
flection of a quasi-particle rather than a free electron, 
because quasi-particles or wave packets composed of a 
complicated combination of plane waves due to U the 
electron-electron interaction manifest themselves in low- 
energy physics. In the one-dimensional Luttinger liquid 
(1DLL), however, the situation is not so simple 0; we 
should ask even the very existence of the Bragg reflection 
and this constitutes one of the aims of this paper. 

Basically there exist two complementary approaches to 
the many-electron system in a crystal described by the 
Hamiltonian H composed of T the kinetic energy, V, and 
U. One is "the band approach" or (T + V) + U in which 
the problem is reduced to the self-consistent determina- 
tion of an effective one-body potential V by combining V 
with the effect of U in the Hartree-Fock-like mean-field 
approximation. If desired, the correlation effect due to U 
(which is missed in V) can be included by perturbation 
in U with respect to the unpertubed part T+V. Another 
is "the correlated-electron approach" or (T+U) + V in 
which we consider a correlated-electron state defined by 
T+U first and then include V perturbatively. Note that 
these two approaches do not always provide the same 
conclusion, as the discussion on the Mott transition Q 
indicates. We can even imagine situations in which the 
competition between V and U brings about a state which 
neither approach describes well. 

In this paper, an example of such intriguing situations 



is shown on the basis of our finding that the Bragg re- 
flection is usually incompatible with the concept of Lut- 
tinger liquids. More specifically we shall treat the ID 
Hubbard model with an alternating periodic potential at 
half-filling, a system attracting much attention in rela- 
tion to the neutral-ionic transition ||, the ferroelectric 
perovskites [^,^), and the crossover from the band insu- 
lator (BI) at V^>U to the Mott insulator (MI) at V<[/ 
. We have made a study in the density-matrix renor- 
malization group (DMRG) || to obtain the charge and 
spin gaps, A c and A s , as well as the electron localization 
length A with the system size up to 400 sites (which is 
larger than any previous calculations by an order of mag- 
nitude). By comparing these exact results with approxi- 
mate ones given in both the band and correlated-electron 
approaches and also by referring to the very recent result 
of Fabrizio et al. we suggest the appearance of a 

Luttinger-liquid-type quasi-metallic state at the BI-MI 
crossover. Here by "quasi-metallic" we mean that it is 
not the same as "metallic" because of nonzero A c , but 
the state is distinct from either BI or MI by such features 
as very small A c and long A. 

For better illustration of the Bragg reflection of a quasi- 
particle, let us start with 3DEG in a weak periodic po- 
tential V for which H is given in second quantization 
with the plane-wave basis as UM 



H = T+V+U = J2 £kc+ c kCT + E E V(K)cLc k+Ko 

K^O kcr 
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k+qcr C k'-qcr' C k'cr' Ckcr , 



(1) 



q^O kcr k'er' 



where Ckcr annihilates an electron specified by momen- 
tum k and spin a, £k = k 2 /2m-/j with m the mass of 
a free electron and /i the chemical potential, V(K) the 
local electron-ion pseudopotential, and £/(q)=47re 2 /q 2 . 

Single-electron properties can be analyzed by the study 
of the thermal Green's function Gk,k+K('w n ) with lu u a 
fermion Matsubara frequency, defined conventionally as 



l/T 



Gk,k+K(iw n ) = - / ^ T (TVckcr(r)c k+Ko 



(2) 
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In the correlated-electron approach, we consider 3DEG 
without V in the first step. Here Gk,k+K(*w„) is not 
zero only for K = 0. Thus we simply write (iu> n ) and 
this function can be determined by the formally exact 
Dyson equation as shown diagrammatically in Fig. 1(a) 
with use of the vertex function A(k'zuv, \duj n ) defined 
in Fig. 1(c). In the second step, we include V in its low- 
est order by solving the equation in Fig. 1(b) to obtain 
Gk,k+K(*w„) in terms of G^ (iu> n ) and A(k'£uv, kiu; n ). 

In order to make the physics as clear as possible, we 
shall be concerned only with the most interesting case in 
which |K| is equal to 2kp with kp the Fermi wave number 
of 3DEG. Then, a band gap A opens at the Fermi level 
with the value of 2V(K) in the noninteracting electron 
gas, while in the interacting system, by expanding the 
self-energy Y.^ G (iuj n ) [= iuj n — e^ — G^ G (iw„) _1 ] with re- 
spect to Lo n up to first order in line with weak V, we find 
easily that A is given exactly as 2z^ F A(fcpO, — kpO)V(K) 
H with G- kF M F (iLu n ) = z kF A/2[(iuj n ) 2 -{A/2) 2 } where 
Zfc F is the quasi-particle renormalization factor at the 
Fermi surface. Thus the ratio, Zk F A(kpO, —kpO), quanti- 
fies the effect of correlation on the Bragg reflection. 

We estimate this ratio quantitatively by employing 
the local-field correction G+(q) to represent the effect 
of the vertex function. Using the values of G+(q) as 
supplied by quantum Monte Carlo ]l3]| , we can calcu- 
late the ratio, together with its components, z kF and 
A(fc F 0, -fcpO) = {l + ar s [l-G + {2k F )]/2n}- 1 , as a func- 
tion of r s the interelectron distance in units of the Bohr 
radius with a = (4/97r) 1 / 3 «0.521. The results are shown 
in Fig. 1(d) in which we see that A(fcpO, — kpO) sw 1. 



(a) In the Electron Gas 



G k G l. 
(b) In a Weak Periodic Potential 

X V(K'-K) 



(d) Band Gap 




FIG. 1. Diagrammatic representation of Gk^+n(iui n ) (a) 
without V and (b) with it in its lowest order. The vertex 
function A(k'iw n / , kiw n ) is defined in (c) with use of the 
irreducible electron-hole effective interaction I. In (d), ef- 
fect of correlation on the Bragg reflection as quantified by 
A/2V(K) = z feF /{l + Qr s [l-G+(2fc F )]/27r} is shown. 

Deviation of the vertex function from unity accounts 
for the many-body effects on the electron-ion interaction 
Vei— ion and physically this should be small at metallic 
densities (l<r s <6) for q — 2kp corresponding to the in- 
terelectron distance in real space; V^i-ion is not modified 
much from the bare one in the neighborhood of an elec- 
tron at this distance due to the absence of other electrons 
by correlation. Thus the ratio is essentially determined 
only by z kF or the weight of the coherent part, confirm- 
ing a naively anticipated result that the Bragg reflection 



occurs only in the coherent part of a quasi-particle. 

The importance of z kF in the Bragg reflection prompts 
us to investigate 1DLL in which Zk F vanishes. The Hamil- 
tonian is basically the ID version of Eq. ([[]); is lin- 
earized as Ek — Vpfo — kp) [i>f(— k — k-p)] for the right- 
[left-]moving branch in T and coupling constants, g±, g 2 , 
and (74, corresponding to the backward scattering, the 
forward scattering between the opposite branches, and 
the forward one within the same branch, respectively, 
are introduced in U ]l5| ]. Using a ka [bka] the annihila- 
tion operator for an electron in the right- [left-] moving 
branch, we can write V as 



V = vJ2(a 



2k F <7 



+ 6 



O-kcr), 



(3) 



in which only the K — 2kp part for V(K) is retained. 

The assertion we shall make is that A c due to V cor- 
responding to A in the Fermi liquid vanishes in 1DLL. 
Our strategy to prove it is to find a criterion as to when 
V turns out to be an irrelevant perturbation. Since the 
key quantity is (a£ Fa b- kF(T ) or TJ2 U „ G-k F ,k F ( iuj n), we 
evaluate the expectation value by treating V as a linear 
perturbation to the system decribed by the Hamiltonian 
T+U. Then, the Kubo's formula provides us 



(H F .b-k F a) = ■= limiV(2fc F ,^) 



(4) 



where N(q,u>) is the charge-density response function. 
Thus the problem is reduced to evaluating N(2kp,u>) in 
the uj~^0 limit for the system of T+U 14 1. The behav- 
ior of this function is known well jy| and the result is 
N(2kp,oj — > 0) oc lo 1 "' 1 with 7 P , given by 



7 P ; 




1 + (.94 + .9i - 2g 2 )/2iTVp 



9i 



2g 2 )l2-KVp' 



(5) 



for gi > 0. This leads us to conclude that (a^ Fa b^k F a) 
vanishes for gi > 2gi even in the presence of V. More gen- 
erally, the expectation value is zero in the phases charac- 
terized by N(2kp,uj — > 0) = 0, indicating the irrelevance 
of V . This irrelevance implies that A c due to V vanishes, 
because the system is the same as that with v = 0. 

Physically low-lying excitations in 1DLL are rigor- 
ously represented by sound waves with wavelengths much 
longer than l/2kp. Thus the effect of V manifests itself 
after the average of V over a distance longer than its pe- 
riodicity, which is null and leads to the complete absence 
of the Bragg reflection. This statement ceases to be valid 
if V is large enough to destroy the Luttinger-liquid state 
itself. In this sense, the concept of the Bragg reflection 
is incompatible with that of Luttinger liquids. 

So far no lattice periodicity is considered in T+U and 
thus the concept of electron filling is irrelevant. Now we 
include it by treating a ID system at half-filling on the 
lattice prescribed by T+U with V possessing periodicity 
of two lattice units. In site representation, H is given by 



2 



H = T+V+U= -t (4„c j+la +c+ +la c j(T ) 
jo 

+ v Yli.- 1 ) 3 eta !* + u Yl c n c n c ti c Ji- ( 6 ) 

JCT j 

For the study of competition between V and U in the 
whole region from w>u to v <Cu, we implement DMRG 
to calculate E(N^,Ni) the ground-state energy with N a 
the fixed number of cr-spin electrons in the L-site sys- 
tem under the open-boundary condition. At size L, the 
charge and spin gaps, A C (L) and A S (L), are given as 



A e (L)=£7(| + l,|)+2^-l,§)- 
A S (L)=E(- + 1,--1)-E(-,-). 



(7) 

(8) 



By using a finite-size scaling as Aj(Z>) = (A 2 +Ai/L 2 + 
/ L 3 + • ■ ■) x / 2 for i = c or s M , we extrapolate the data at 
L = 50, 100, 200, and 400 to obtain the values at L = oo, 
A c and A s . In Fig. 2, the results thus obtained are plot- 
ted as a function of u at v = 0.5t. As u increases, both 
A c and A s decrease from 2v (the band gap at u = 0) and 
become very small at around u = 2.6t. With the further 
increase of u, A c increases very rapidly, while A s remains 
to be zero. In the inset of Fig. 2, the gaps near u = 2.6t 
are shown in detail. The same overall behavior of the 
gaps is seen for other values of v of the order of t. 




3 4 5 
u (units of t) 

FIG. 2. Charge and spin gaps as a function of u at v = 0.5t. 
Approximate values in the band and correlated-electron ap- 
proaches are shown by dashed and dotted curves, respectively, 
while those in HF by the dotted-dashed curve. The inset 
shows the results around u — 2.6t in a magnified scale. 

Let us analyze these results by the comparison with 
those in both band and correlated-electron approaches. 
The former approach begins with the Hartree-Fock (HF) 
approximation which amounts to the^one-body problem 
described by T+V where we define V in the same form 
of V in Eq. (^) by replacing v into v, determined through 



v = v 



* dk 



-2L 2n Vw 2 +4£ 2 cos 2 k 



(9) 



In HF we obtain the band insulator (BI) brought about 
by the full Bragg reflection with A c = A s — 2v. Since v 
is always positive, these gaps never vanish, which clearly 



contradicts the exact results for large u. Even for u as 
small as t, they do not agree well with the exact ones, as 
seen by the dotted-dashed curve in Fig. 2. 

However, we achieve a surprisingly good improvement 
by including the correlation effect in second-order per- 
turbation in U with the unperturbed basis in T + V, as 
shown by the dashed curve in the inset of Fig. 2. In fact, 
the exact gaps are reproduced quite accurately in this 
(T+V) +U approach for u from up to about u c o = 2Aht. 
Note that the correlation effect included in this way does 
not separate A c from A s . Physical mechanism to reduce 
the gaps from those in HF is the same as explained in 
3DEG, namely, the reduction of the the Bragg reflection 
on the conversion from a free electron to a quasi-particle. 
Thus we conclude that the state realized in the system 
for u<u c o is the correlated BI. 

For u larger than about w C 2 = 2.90i, on the other hand, 
A s vanishes, while A c does not, indicating that the state 
is the MI. The effect of v can be included in A c in the 
(T+U) +V approach; by examining the exact solution in 
T+U |p^| , we can deduce an expansion for A c in u^ 1 up 
to third order as 

. f 2 y 2 -\-4-v 2 f 4 

A c « u -2vV+4t 2 + 81n2 — 6C(3)^, (10) 

with C(3) w 1.202. The above result is plotted by the dot- 
ted curve in Fig. 2 in which we see that the exact result 
is reproduced quite well for u larger than 5i. 

Now we need to clarify the nature of the state for u at 
the BI-MI crossover, ranging from u c q to u C 2- For that 
purpose, we calculate Dl & dimensionless localization pa- 
rameter introduced by Resta and Sorella [0 as 



D, =-Lhx 



(11) 



with ^ the ground-state wavefunction at iVj =N± =L/2 
under the open-boundary condition for the L-site system 
and Xj the position operator at site j. Extrapolation of 
Dl to the L — ► oo limit gives the value D which is related 
to A through A = \[D /2tt in units of the lattice spacing. 

The calculated results for both D and Dl at various 
L's are shown as a function of u in Fig. 3(a). For either 
u < u c q or u > u C 2, we see that Dl converges to D at L 
as small as 100, while for u in-between, even L = 400 is 
not large enough for the convergence. In particular, D 
seems to diverge for u around u c i = 2.65£. Divergence in 
D, or equivalently that in A, implies the appearance of a 
metallic state, indicating the vanishment of A c . 

As for A c — at u = u c x and the state for u c i < u < u C 2, 
Fabrizio et at. ]l0[ ] have suggested a spontaneously dimer- 
ized insulating phase (SDI) by analytical arguments. In 
fact we find that A decreases quite rapidly to be less than 
twice the lattice spacing as u increases from u c x, implying 
an insulating behavior. Thus we conclude that an insu- 
lating phase, most likely SDI, is realized for u c \ <u< u C 2- 
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(b) Size Dependence 
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(a) Dimensionless localization parameter Dl as a 
function of u at v = 0.5£ for various system sizes, together 
with the value extrapolated to L = oo. The values of Dl 
corresponding to A = l, 2, and 3 are also indicated, (b) Size 
dependnece of both A C (L) and A S (L) at u = 2t, 2.6i, and 3t. 

A remaining problem is that A c never becomes zero in 
Fig. 2, although we now know that it should be zero at 
least at u = u c \. Thus we reexamine the size dependence 
of both A C (L) and A S (L) carefully and the typical re- 
sults are plotted in Fig. 3(b). Let us first analyze them in 
terms of A c — A s or the spin-charge separation. Because 
of one dimensionality, one might assume that A c =f= A s 
should be the case as long as it^O, but this is not true; 
at u = 2t in the BI region, we find that both A c and A s 
coincide up to at least five digits (which exceeds numeri- 
cal accuracy) at L = 200 or larger. On the other hand, we 
know that A c — A s = — A s ^0 at u — u c \. Therefore there 
definitely exists a value of u at which A c — A s begins to 
deviate from zero. We identify u c0 as such a value and 
thus only for u larger than u c o the IDLL-like spin-charge 
separation occurs. In this sense, we consider that it c n 
gives a sharp phase boundary. 

In order to find more detailed features of the state at 
u from u c o to u&, let us look at Fig. 3 again. As rep- 
resented at u = 3t in Fig. 3(b), A S (L) converges to zero 
very nicely with the increase of L in MI. Similar conver- 
gence of both charge and spin gaps is obtained in SDI as 
well. However, a distinct behavior is seen in the gaps for 
u c o < u < u c \ as illustrated at u — 2.6t; both A C (L) and 
A s (L) at L = 400 seem to be larger than those extrapo- 
lated from the data at smaller L, implying that our data 
for the gaps are not accurate enough for these it's. This 
inaccuracy should be due to the large A which is always 
longer than the two lattice units for these it's as indicated 
in Fig. 3(a). This difficulty in obtaining exact values for 
the gaps in this region can be overcome only by a more 
accurate calculation of energies at much larger L. Such 
a calculation is not feasible at present, but we can safely 
conclude even at the present time that A c is very small, 
i.e., much smaller than O.lt at most of the values of u in 
this phase. The nature of small A c and long A suggests 
us that a Luttinger-liquid-type quasi-metallic phase ap- 



pears for u c q < u < u c \. Incidentally, each electron in this 
phase feels V with the spatial average over A which is 
longer than the periodicity of V . This implies that the 
effect of V on the electrons is very small and thus this 
should be the reason why the feature of BI is lost in the 
state for u c q < u < u c x- Here again we find that the Bragg 
reflection, a crucial concept to define BI, is incompatible 
with the Luttinger-liquid feature. 

Finally we note that both this phase and SDI deserve 
special attention, because they are not anticipated in 
both band and correlated-electron approaches; their ex- 
istence is entirely due to the competition of V and U. 

In conclusion, we have investigated the effect of elec- 
tron correlation on the Bragg reflection in a variety of 
situations in various approaches. We have found the in- 
compatibility of the Bragg reflection with the Luttinger 
liquid, based on which we have suggested a Luttinger- 
liquid-type quasi-metallic state at the crossover from BI 
to MI via SDI. 
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